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The longitudinal double spin asymmetry All in the reaction pp ttX has been measured at 
00 ' RHIC with extremely interesting consequences. If the gluon polarization in a proton were as big as 

I needed to resolve the famous "spin crisis" then All would be large and positive. Latest RHIC results 

>— ^ , indicate that ^ll is small and disfavour large positive values of the gluon polarization. We examine 

whether the soft mechanisms (Collins, Sivers, Boer-Mulders) , essential for generating transverse 
^ ' single spin asymmetries, have any significant influence on All, and whether they could alter the 

Q , conclusion that the gluon polarization is necessarily small. It turns out that the contribution from 

■ these effects is essentially negligible. 

O ' PACS numbers: 13.88.+e, 13.85.-t, 13.85. Ni 
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I. INTRODUCTION 



Large transverse single spin asymmetries (up to 40%) have been observed in a multitude of reactions for over three 
Qj^' decades, whereas such asymmetries are tiny (< 1%) in the standard leading twist QCD parton model. To explain 
] the size of these asymmetries Sivers and Collins [1, 2] introduced new soft mechanisms, utilizing, as an essential 
ingredient, the intrinsic transverse momentum of partons [3]. Later other similar mechanisms were shown to be 
I possible, for example the Boer-Mulders mechanism [4, 5]. 
^ . Although these mechanisms were invented in order to produce transverse asymmetries, it turns out that they also 
^\ ' contribute to the longitudinal double spin asymmetries and to the total cross-section [6]. For the latter, it has been 
\^ , shown that the effect of the soft functions is negligible. However, it was found that intrinsic transverse momentum 
•/^ ' per se significantly affects the value of the cross-section [7]. That this should be the case at lower energies had already 
'""J \ been noticed by Field and Feynman [8] , and later by Vogelsang and Weber [9] on the grounds that taking into account 

■ intrinsic fcj^ is a particular method of including higher twist corrections. 
^-H For the longitudinal asymmetries, the question is much more delicate, for the following reason. One of the most 

. important reactions measured at RHIC is the double spin longitudinal asymmetry All, which has been found to be 
very small, and which has been used, based on a leading twist coUinear treatment, to confirm the growing belief that 
the gluon polarization is far too small to explain the "spin crisis in the parton model" [10]. 

Given that All is so small, and that the implications of this are so important, we felt it necessary to check whether 
the soft mechanisms can have a significant impact, in particular whether they could influence the above conclusion 



' about the gluon polarization. We have also checked for any significant sensitivity in All to intrinsic transverse 
momentum. 

The plan of the paper is the following. In Section II we briefly recall the formalism used for our calculation, which 
includes the full non-collinear kinematics of the scattering process. In Section HI we present the "kernels" for the 
calculation of each partonic contribution to the polarized cross-sections. In Section IV we show and discuss our 
phenomenological results for the longitudinal double spin asymmetry in inclusive neutral pion production at RHIC. 
Finally, in Section V we draw our conclusions. 



II. FORMALISM 



Here we simply sketch the main aspects of the formalism; for details of the approach we refer to [6] . The longitudinal 
double spin asymmetry All for the reaction pp ttX is defined as 

_ da++ - da+- _ da++ - da+- 

~ d<T++ + da+- ~ 2d&^P ' ^ ' 

where the labels refer to the helicities of the protons. 
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The general expression for the differential cross-sections for the polarized hadronic process {A, Sa) + {B, Sb) C+X 
is given by 

£;pd^(A,S.4) + (B,SB)-C+X r dXadXbdz 



d^Pc 



a,fc,c,rf, { A} 

X Ma,,a,;A„,a, Ma1,a,;A'„,a; + i + ^) Dl^::xnz, k^c) , (2) 



which involves a (factorized) convolution of all possible hard elementary QCD processes, ab — > cd, with soft partonic 
polarized distribution and fragmentation functions. In Eq. (2) s, t and u are the Mandelstam variables for the partonic 
reactions. The detailed connection between the hadronic and the partonic kinematical variables is given in full in 
Appendix A of Ref. [6]. A discussion of some technical details, like, e.g., the infrared regulators related to small 
partonic scattering angles, can be found, for example, in Ref. [7]. 

Let us simply recall here, for a better understanding, the physical meaning of the different factors in Eq. (2): 

• pll^y''^ is the helicity density matrix of parton a inside the polarized hadron A, with spin state Sa\ it describes 

the parton polarization. fa/A,SAi^a, k±a) is the number density (or distribution) of unpolarized partons a inside 
the polarized hadron A, Sa- each parton carries a light-cone momentum fraction Xa and a transverse momentum 
k±a- Similarly for parton b inside hadron B with spin Sb- 

• The polarized cross-sections for the elementary partonic process (a, Sa) + (b, si,) — > (c, Sc) + d are expressed in 
terms of products of the helicity amplitudes Ma^.a^;A^.Aj^- 

• The factor D^'^y'^ (z, k±c) describes, again in the hehcity basis, the fragmentation process c C+X, according 
to which a polarized parton c fragments into an unpolarized hadron C carrying a light-cone momentum fraction 
z and a transverse momentum kj^c- 

• J{k±c) is a kinematical factor, numerically very close to 1 for RHIC kinematics. All details can be found 
in Ref. [6]. Throughout the paper, we work in the AB cm. frame, assuming that hadron A moves along the 
positive Zcm-B^xis and hadron C is produced in the (XZ)cm plane, with {pc)x^rn > 0- 

Eq. (2) is written in a factorized form, separating the soft, long distance from the hard, short distance contributions. 
The hard part is computable in perturbative QCD, while information on the soft one has to be extracted from other 
experiments or modeled. As already mentioned and discussed in Refs. [6, 7], such a factorization with non-coUinear 
kinematics has never been formally proven. Indeed, studies of factorization [11-14], comparing semi- inclusive deep 
inelastic scattering (SIDIS) and Drell-Yan reactions have indicated unexpected modifications of simple factorization, 
and the situation for single inclusive particle production in hadron-hadron collisions is not yet resolved. Thus, our 
approach can only be considered as the natural extension of the coUinear case and a reasonable phenomenological 
model. Of course, the perturbative calculation of the hard part is only reliable if the hard scale - in this case the 
square of the transverse momentum of the final hadron, - is large enough. It turns out that the data on unpolarized 
cross-sections in hadronic collisions at low-intermediate energy scales suggest [7] an average value of fc^. = I^J-P — 0-64 
(GeV/c)^ for the intrinsic transverse momentum of the parton distributions. On the other hand, both unpolarized 
hadronic cross-sections at RHIC energies [15] and the Cahn effect in SIDIS [16] are rather well reproduced by using 
(fc^) ~ 0.25 (GeV/c)^. We shall therefore study how the contributions to All depend on the value of (fc^^). 



III. KERNELS 



As we can see from Eq. (2), the computation of the cross-section corresponding to any polarized hadronic process 
{A, Sa) + {B, Sb) C + X requires the evaluation and integration, for each elementary process a + b ^ c + d, of the 
general kernel 

^(^A, S's)"''^'"^ = ^Pxi'^xf'^ fa/A,SAi.^a,k±a)Px[^X^'' fb/B,SBiXb,k±t) 
{A} 

xAfA,,A,;A„.A,A>A^ ,A,;A^,a; ^t^''^^' '^^^^ ' 
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While the hadronic process {A, Sa) + {B, Sb) C + X takes place, according to our choice, in the {XZ)cm plane, 
all the elementary processes involved, A{B) a{b) + X, ab cd and c ^ C + X do not, since all parton and hadron 
momenta, p^, p^, have transverse components k±a, fe±h, k±c- This "out of (XZ)cm plane" geometry induces 
phases in the fragmentation process, in the distribution functions and in the elementary interactions, which have to 
be taken into account. Thus, the independent helicity amplitudes for the elementary pQCD processes ah — > cd, with 
massless partons, can be written as [6] 



+ ,+;+, 



-;+,- 



where the amplitudes Mi 2 3 are the real planar amplitudes defined in the partonic ab cd cm. frame, 
A/i" = Ml ^ = a7" _ Ml = All _L = All _ All = All ^.^ _ = Ml 



,+;+■ 



-,+;-,+ 



'-,+;+.- 



(4) 



(5) 



as required by parity invariance. The phases </'i,2,3 are complicated functions of the polar and azimuthal angles of the 
transverse momenta, fc_LQ, fc^h and fc^c"? and their explicit expressions can be foimd in Rcf. [6]. The relations 



+ ,+;+,+ 



A/* 



A/* 



,+;+, 



(6) 



follow from Eqs. (4), (5) and from the fact that the phases tfi change sign by helicity inversion [6]. Note that the + 
and — subscripts refer to (+1/2) and (—1/2) helicities for quarks, and to (+1) and (—1) helicities for gluons. There 
are eight elementary contributions ab — > cd which we have to consider separately 



qaqb QcQd , gagt gc9d ■. 
qg-^ qg, gq gq , 
qg^ gq, gq^ qg, 

gagb qq, qq-^ gcgd , 



(7) 



where q can in general be either a quark or an antiquark. The subscripts a, b, c, d for quarks, when necessary, 
identify the flavour (only in processes where different flavours can be present); for gluons, these labels identify the 
corresponding hadron (a — > A, b —f B, c C). By performing the explicit sums in Eq. (3), we obtain the kernels for 
each of the elementary processes. Note that the new aspect of our calculation is the appearance of the phases which 
is a reflection of the non-coUinear kinematics. 

The computation of the denominator/numerator of All in Eq. (1) requires the evaluation of the kernels [S(+, +) ± 
S(+, — )] respectively. The expressions for the sums of kernels, which are relevant for the unpolarized cross-section, 
are given in Ref. [6]. Here we give in detail the expressions for the differences. They are calculated from the general 
kernel given in Eq. (3). In the following certain terms are underlined: these are terms which vanish after integration 
over the angles of the momenta k±a, l^±b, fc±c in Eq. (2); we shall further comment on that at the end of this Section. 
(f>Q is the azimuthal angle of the hadron C in the parton c helicity frame and its expression in terms of the angles of 
k±c is given in Appendix A of Ref. [6] . Notice that all angular dependences of the kernels are explicitly extracted and 
the parton distribution (PDF) and fragmentation (EE) functions only depend on the magnitudes of the transverse 
momentum vectors. 

• qaqt ^ qcqd contribution 

[!](+,+) -S(+,-)]9='?^^«=«^ = 

Afl/^ixa,k^a)Afl^4xb,k^b) \ \M?\' - \Ml\' ~ |Af°pj Dc/ciz,k^c) 



/^(Xa, kj_a) Af^,^{xb, kj_b) C0S((P3 - (^2) 



+AfsjAixa, k±a) Aflj_^{xb, k±b) sin(v33 - 1^2) 



i2MlMl) Dc/c{z,k^c) 



fa/AiXa, fc±a) Af^^^ixb, k^b) M? M^ sin(^i -^3+ (/)g) A^^c7/,T (z, fc±c) ■ 



(8) 



Notice that we have used the relations Af^^^^{xa, k±a) = Af^^^^{xa, k^a) and fa/+{xa,k±a) = fa/Aixa, kj_a), 
see Appendix B of Ref. [6]. The channels qq — > qq etc. arc formally identical to qq qq with amplitudes defined 
properly in Rcf. [6] . 
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gq — !■ gq contribution 



Mlj+{xa,k^a) Mlj+{^b, k^b) [ - lAfOp] Dc^giz, k^c) ■ (15) 

The physical content of the above expressions is interesting. First note the complete formal symmetry between the 
qq — > qq kernel in Eq. (8) and the gg —^ gg kernel in Eq. (9). These kernels contain the largest variety of contributions, 
and the kernels for all the other partonic processes can be formally read off from these by the suppression of certain 
terms. 

In the second line of both expressions, Eq. (8) and Eq. (9), we recognize the product of the fcj_-dependent helicity 
distributions, A/^ /^{xq, k±q) = Ag(xq, k±q) and A/^ ^^g) = '-^Qi^gt k±g) for quarks and gluons respectively, 

and the unpolarized fragmentation function Dc/c{z, k±c)^ with no azimuthal phases. In the third line of Eq. (8) we 
have two parton distribution functions, A/^ /^{x,k±), referring to quarks transversely polarized, along the x-axis, 
inside longitudinally polarized nucleons, coupled to the unpolarized fragmentation function. Analogously, in the third 
line of Eq. (9) we have two parton distribution functions, Af^^^_^{x, kj_), which are related to the linear polarization 
of a gluon inside a longitudinally polarized nucleon. Correspondingly the fourth line of Eq. (8) refers to one quark 
transversely polarized along the cc-axis inside a longitudinally polarized nucleon and the other, A/J y^(a;, fcj_), trans- 
versely polarized along the y-axis inside an unpolarized nucleon - the latter is the Boer-Mulders function ~ coupled 
to the unpolarized fragmentation function. Analogously, in the fourth line of Eq. (9) we have Af^_^^^{x, k±) and the 

"Boer-Mulders-like" gluon function, Af^^^^(x,k±), referring to a linearly polarized gluon inside an unpolarized nu- 
cleon. For a more complete explanation of the physical meaning of these functions see Appendix B of Ref . [6] . Finally, 
the last line of Eq. (8) contains the Collins fragmentation function, Dq/^t {z, k±c), coupled to an unpolarized 
parton density and a transversely polarized one. In the case of the gluon, in the last line of Eq. (9), there appears a 
giuonic analogue of the Collins fragmentation function, A^ D(j fq-s{z, k±c), describing the fragmentation of a linearly 
polarized gluon into an unpolarized hadron. 

Ignoring the underlined terms which vanish upon integration, we see that compared to the standard coUinear ap- 
proach, we have extra contributions involving quarks polarized transversely along their z-axis in a longitudinally 
polarized nucleon, appearing in Eqs. (8), (10) and contributions involving linearly polarized gluons inside a longi- 
tudinally polarized nucleon, appearing in Eqs. (9), (11). Notice that the processes in Eqs. (12)-(15), initiated by 
quark-gluon elementary scattering, get contributions only from the usual terms, which survive in the collinear case. 

The demonstration of the vanishing upon angular integration of the underlined terms in Eqs. (8)-(15) requires 
a detailed study of the kinematics and of the relationships between the angular integration variables appearing in 
k±a, k±c in Eq- (2) and the phase variables (pi.2,3 and 4>^ [6]. We have also numerically checked that this is 
indeed the case. 

Notice that a parity transformation implies ipi —ipi [i — 1,2,3) and 0^ Thus the odd sin terms in 

Eqs. (8)-(15) must vanish if parity is conserved. 

Another simple, but interesting example of such a vanishing can be obtained by considering, within the same 
formalism, the expression of the kernels for the longitudinal single spin asymmetry A^, which we know must vanish 
in a parity conserving theory. The kernels themselves are not zero, but under integration do vanish. This is another 
very stringent test of the correctness of our formalism. For A^^ for the partonic channel qaqb ~* QcQd, we have for the 
numerator of the longitudinal single spin asymmetry the following expression: 

[S(+,0)-I](-,0)]''"«''^'?'=«<' = 

Af:j4xa,k^a)Af'^js{xb,k^,) {2M°M°) sin(y.3 - ^2) i^c/c(^, fc±c) 

~Afl/^{xa,k^a) h/B{xb. k^b) sin((^i - ^2 + ^g) M'^ A4 A^Dc/c^ {z, k^c) , (16) 
and, again, all terms - being odd functions of (pi and (j)^ - vanish, as they should, upon angular integration. 
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IV. PHENOMENOLOGY: All AT RHIC 

The longitudinal double spin asymmetry All for inclusive neutral pion and jet production in proton-proton scat- 
tering at y'i = 200 GeV has been measured at RHIC in various runs, respectively by the PHENIX [17-19] and 
STAR [20, 21] Collaborations. The first pubhshed PHENIX experimental data [17] showed results for A^^ at mid 
rapidity compatible with negative values. This was quite puzzling, since is a positive quantity in the collinear 
parton model [22], at least at low px where it is dominated by gg gg elementary scattering processes, see Eq. (9) in 
which it can be shown that |M°|^ — [Mj |^ — l-^s [^ > 0. More recent and precise data from both collaborations [18, 20] 
exclude the possibility of a large and negative Al^: in two subsequent RHIC runs, Run 5 [19, 21] and Run 6 (results 
from Run 6 have only been presented as "preliminary" [23, 24]), they confirm and reinforce the statement that All 
is very small and compatible with zero over the whole px range covered. 

An earlier comparison of present RHIC data with collinear next-to-leading order (NLO) QCD calculations of 
All [25] disfavoured large positive values for Ag, definitely excluding scenarios where Ag is as large as the unpolarized 
gluon distribution function, g, at low scale. Instead the data were in better agreement with the predictions obtained 
by assuming Ag = or even Ag ~ —g at the initial scale [26]. A recent statistical analysis shows that the PHENIX 
Run 5 data are compatible with both Ag — and the "standard" GRSV parametrization [27], while it rules out the 
Ag = —g (at the initial scale) hypothesis [19]. A newest update of this analysis, which includes the preliminary data 
from PHENIX Run 6, favours the Ag = scenario over the standard GRSV [24]. (Note that in this Section we have 
adopted the common, short-hand notation A/' = Aq and Af^ = Ag for the helicity distribution functions, 
while fq/p = q and fg/p = g for the unpolarized distribution functions, for quarks and gluons respectively). 

Our goal is to explore whether the new mechanisms permitted by the presence of partonic intrinsic transverse 
momenta, obtained in a general and fully non-coUinear kinematics, could affect the above conclusions, which are 
based on the analysis of All in the collinear configuration, i.e. taking into account only the terms proportional to 
Aq{x) and Ag(x). Could the "new" contributions shown in Eqs. (8)-(ll) turn the longitudinal double spin asymmetry 
All into a very small (or even slightly negative) quantity without the need to assume Ag to be zero or negative? 

We have studied All at RHIC, for the PHENIX kinematics, ^/i = 200 GcV and |?7| < 0.35 (numerical calculations 
are performed at 77 = 0) and evaluated each separate contribution to All, according to Eqs. (8)-(15). Since we have 
no knowledge of the parton densities Af^ and Af^^^_^ we maximized them in order to see whether, in principle, 
they can have a significant effect on All- We thus used for them the corresponding unpolarized parton densities and 
adjusted the signs so that all contributions add up coherently. 

For the hehcity distributions we have used the sets GRSV2000 [27] and LSS05 [28]. The unpolarized cross-section 
and the maximized contributions to the numerator of All have been calculated using the GRV98 set [29] and the 
MRSTOl set [30] respectively. For the fragmentation functions we have used the KKP set [31] and, for comparison, 
the Kretzer set [32]. The transverse momentum dependence has been included by means of a factorized Gaussian 
smearing, for all the parton distribution and fragmentation functions 

f(x,k^) = fix) , (17) 

D{z,k^c) = D{z) , (18) 

with a constant and flavour independent parameter Vv^uT = ^Oj assumed to be the same for all quark flavours and 
for gluons; we shall study the effect of changes in the value of ko. Guided by our previous work, we compared the 
results obtained using three different values for kg: kg = 0.8 GeV/c from studies on the unpolarized pp scattering 
cross-sections and single spin asymmetries [7], fco = 0.5 GeV/c from fitting the Cahn effect in SIDIS [16], and fco = 0.01 
GeV/c to recover the collinear configuration. For the fragmentation functions, we take {k\(j) = {k\) everywhere. We 
have checked that variations in (fc^^) induce negligible changes in All- 

It turns out that the new non-coUinear soft contributions containing the PDFs A/J y_|_(x, A:^) and Af^^^^{x^k±^), 
even if maximized, are totally negligible. In fact, in the RHIC kinematical regime considered their maximized contri- 
bution docs not exceed, in the lowest pT range, few percent of the usual terms (already present in the collinear case), 
becoming much smaller at larger pr- We have checked that this result remains true also at lower energies. Let us 
remark that a similar situation holds also for the unpolarized cross-section [7]. Although the two additional terms 
(with respect to that already present in the collinear case) are of course different in this case [6] , involving respectively 
the convolution of two Boer-Mulders functions with an unpolarized fragmentation function and the convolution of a 
Boer-Mulders function and an unpolarized distribution with the Collins fragmentation function, their total maximized 
contribution reaches at most 1% of the usual term, being even smaller on the average. 
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Figure 1: All for the process pp tt^X at y/s — 200 GeV and rj = 0, pfotted as a function of pr, calculated with different 
choices of {k'j_) = ko in the PDF/FFs, compared to PHENIX data, Run 5 [19]. The solid line corresponds to the choice 
ko = 0.01 GeV/c in both PDFs and FFs. The dashed line corresponds to ko = 0.8 GeV/c in PDF/FFs. The PDF sets are 
LSSQ5 [28] and MRSTOl [30], the FF set is KKP [31] and the factorization scale is Q = pr- Notice that the changes in All 
induced by varying the value of (fcx) much smaller than those obtained by choosing different sets of distribution functions 
and/or factorization scales, see Fig. 3. 

Therefore, we conclude that there is no way for the extra contributions induced by the presence of partonic intrinsic 
transverse momenta to alter the size oi All- We have checked that this conclusion is not sensitive to the choice of the 
mean intrinsic transverse momentum ko. In fact, Fig. 1 shows that in general All depends very little on the different 
choices of ko; in particular, All decreases when increasing the width of the gaussian, but compared to data this 
variation is quite negligible. This result can be understood because the k± dependence is given by the same gaussian 
for all distribution and fragmentation functions and at mid rapidity the M amplitudes depend very mildly on k±. 

It is interesting to notice that the corresponding unpolarized cross-section is also almost independent of the value 
assigned to the average intrinsic transverse momentum fco , while it turns out to be more sensitive to the choice of the 
factorization scale and of the fragmentation function set, as we show in Fig. 2, consistently with the NLO collinear 
pQCD calculations. The comparison with PHENIX data [19] is well satisfactory. The solid lines correspond to the 
factorization scale Q = px, the dashed lines to Q = pt/2, using the GRV98 [29] PDF set and the KKP [31] (thick 
lines) or the Kretzer [32] (thin lines) FF sets. Results for the STAR and BRAHMS kinematical regimes at y/s = 200 
GeV can be found in Rcf. [15] and show similar agreement with data when adopting the same average fcj^'s as in 
Ref. [16]. 

Contrary to what happens for the k± dependence. All is sensitive to the choice of the PDF set and of the scale. 
In Fig. 3 we show All calculated in an almost collinear configuration, a/ (fc^) = 0.01 GeV/c, and for two choices 
of scale, Q = Pt and Q = pt/"^- Using the LSS05/MRST01 PDFs, the dotted line corresponds to Q = Pt/2 and 
the dash-dotted line to Q = pr- Using the GRSV2000/GRV98 PDFs, the sohd line corresponds to Q = pt/2, and 
the dashed line to Q = pr- As can be seen, the variations induced by different choices of PDF sets and scale are 
quite large, larger than those produced by changes in the fco value; nevertheless, all these curves are compatible with 
present experimental data (we have checked that these same conclusions hold also when adopting Q = ^pt). However, 
very precise data on All in the future might be able to distinguish between various sets of PDFs. Data collected at 
different energies [23] will also be very useful to cover presently unexplored regions of the Bjorken x variable. 

Concerning the dependence of All on the set of fragmentation functions, we have checked, adopting again the KKP 
and Kretzer sets, that this is almost negligible over the whole pr range considered. Only at the largest px values, 
where All data show large experimental errors, there is some residual dependence. This result can be understood, 
since both the numerator and the denominator of All contain the unpolarized FF. 

V. CONCLUSIONS 

We have examined, at leading order in perturbative QCD, the effect on the longitudinal double spin asymmetry All 
of allowing the partons to have non-zero intrinsic transverse momentum, and of including in All the contributions 



Vs=200 GeV 
ri=0 
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Figure 2: The invariant unpolarized cross-section for the process pp — > tt^X at ^/s — 200 GeV and 77 = 0, plotted as a function 
of pt, calculated with different FF sets and factorization scales. The thick, solid and dashed lines correspond to the choice of the 
KKP FF set [31], at the factorization scale Q = pr and Q = pt/2 respectively. The thin, solid and dashed lines correspond to 
the choice of the Kretzer FF set [32] , at the factorization scale Q = pr and Q = pt/2 respectively. The PDF set is GRV98 [29] . 
^/(k\) = fco = 0.5 GeV/c for both PDFs and FFs. The experimental data are from the PHENIX collaboration at RHIC, Run 
5 [19]. 




Figure 3: All for the process pp — > 7r°X at y/s — 200 GeV and 7; = 0, plotted as a function of pr, calculated with different 
PDF sets and factorization scales. The solid and dashed lines correspond to the choice of PDFs GRSV2000/GRV98, at the 
factorization scale Q — Pt/2 and Q = pr respectively. The dotted and dash-dotted lines correspond the choice of PDFs 
LSS05/MRST01, at the factorization scale Q = pt/2 and Q = pr respectively. The FF set is KKP. ^(fc^ ) = 0.01 GeV/c for 
both PDF/FF. The experimental data are from the PHENIX coUaboration at RHIC, Run 5 [19]. 



arising from the new soft functions that play a crucial role in transverse single spin asymmetries. The study was 
carried out in the hope that such effects might negate the conclusion that the very small measured values of A^l 
automatically imply that the polarized gluon density is very small. Our analysis indicates that the contribution from 
these effects is negligible and we are forced, at the present stage, to accept the conclusion that the polarized gluon 
density is much too small to explain the "spin crisis in the parton model" [10]. 
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